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Abstract. In this paper a method for embedding cocyclic submatrices with "large" determinants 
of orders It in certain cocyclic Hadamard matrices of orders At is described (t an odd integer). If these 
determinants attain the largest possible value, we are embedding D-optimal designs. Applications 
to the pivot values that appear when Gaussian Elimination with complete pivoting is performed on 
these cocyclic Hadamard matrices are studied. 
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1. Introduction. A Hadamard matrix H of order n is an n x n matrix with 
elements ±1 and HH T = nl. A Hadamard matrix is said to be normalized if it 
has its first row and column all I's. We can always normalize a Hadamard matrix 
by multiplying rows and columns by —1. These matrices must have order 1, 2 or a 
multiple of 4. It is conjectured that Hadamard matrices exist for every n = (mod 4). 
Although no proof of this fact is known, there is much evidence about its validity (see 
[IB] and the references there cited). 

Two Hadamard matrices Hi and Hi are called equivalent (or Hadamard equiv- 
alent, or ii-equivalcnt) if one can be obtained from the other by a sequence of row 
and/or column interchanges and row and/or column negations. The question of clas- 
sifying Hadamard matrices of order n > 32 remains unanswered and only partial 
results are known [15] . 

Problems involving Hadamard matrices sound very easy, but they are notoriously 
difficult to solve. One interesting open problem, among others, is the question of 
the largest pivot encountered during the process of Gaussian Elimination (GE) with 
complete pivoting for an n x n Hadamard matrix H (the so called "growth factor" for 
H). Traditionally, backward error analysis for Gaussian Elimination (GE), see e.g. 
[7J, on a matrix A = is expressed in terms of the growth factor 
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i\j , k = 1,2, ... ,n, that occur during the elimination 
for a choice of pivoting strategy given. Matrices with the property where no row and 
column exchanges are needed during GE with complete pivoting are called completely 
pivoted (CP) or feasible. In other words, at each step of the elimination the element 
of largest magnitude (the "pivot" , denoted by pk) is located at the top left position of 
every appearing submatrix during the process. If A(k) denotes the absolute value of 
k x k principal minor of A, then mathematically A being CP means (or is equivalent 
to) that for each fc, we have that A(k) is greater than or equal to the absolute value 
of any other k x k determinant that includes the first k — 1 rows and columns. This 
is not necessarily the maximum k x k minor of A, but only the maximum kxk minor 
of A when its first k — 1 rows and columns are fixed. 

For a CP matrix A we have 

max{pi,p 2 , • • ■ ,Pn} 



g(n,A) = 



If a matrix is not initially CP, by applying row and column operations with complete 
pivoting we can always bring it in CP form. 

The following lemma gives a useful relation between pivots and minors. 

Lemma 1.1. Jd$ Let A be a CP matrix. The magnitude of the pivots which appear 
after application of GE operations to A is given by 

PJ= A P J \y i = l,2,...,n, A(0) = 1. 



In 1969, Cryer [6] conjectured that if A is a real nxn matrix such that \dij\ < 1, 
then g(n, A) < n, with equality iff A is a Hadamard matrix. In 1991 Gould [14] proved 
that the first part of the conjecture is not true. He found matrices with growth bigger 
than their orders. Thus, the following remains open: 

Conjecture(Cryer) The growth of a Hadamard matrix is its order. 

This conjecture has only been proven for n = 4, 8, 12 and 16 (see [6][9j[22]). Great 
difficulty arises in the study of this problem because TJ-equivalence operations do not 
preserve pivots, i.e. the _ff-equivalent matrices do not necessarily have the same pivot 
pattern. For instance, for n = 16 there are 34 pivot patterns although there are 
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only 5 equivalence classes of Hadamard matrices for this order. Furthermore, many 
pivot patterns can be observed by permuting the rows and columns of any 20 by 20 
Hadamard matrix and there are just 3 inequivalent matrices. 

However, the existence of D-optimal designs (and other specific submatrices with 
concrete determinants) that exist embedded in a Hadamard matrix have provided 
some clues on the pivot patterns (see [26 ] I24 j ). 

A D-optimal design of order n is an n x n (1, — l)-matrix having maximal deter- 
minant. Here and throughout this paper, for convenience, whenever a determinant or 
minor is mentioned, we mean its absolute value. The question of finding the deter- 
minant of a D-optimal design of order n is an old one which remains unanswered in 
general. 

In 1893 Hadamard proved in [15] that for every (—1, l)-matrix M, 
(1.1) det(M)<n%. 



We recall that the original interest in Hadamard matrices stemmed from the fact 
that these matrices are the only ones that satisfy equality in ([1,1)1 , 

This has led to further study and tighter bounds for the maximal determinant for 
all (-1, l)-matrices of order n ^ (mod 4) have been found (see EH [HI [28] [20]). 
For instance, when n = 2 (mod 4), Ehlich in [10] and independently Wojtas in [28] 
proved that 

(1.2) det(M) < (2n-2)(n-2)T^. 



In order for equality to hold, it is required that there exists a (— 1, l)-matrix M of 

' where L = (n — 2)J« + 2 . Here, as usual, /„ 



order n such that MM T = 



L 

denotes the identity matrix of order n, and J n denotes the n x n matrix all of whose 
entries are equal to one. In these circumstances, it may be proven that, in addition, 
2n — 2 is the sum of two squares, a condition which is believed to be sufficient (order 
138 is the lowest for which the question has not been settled yet, [13]). 

In the early 90s, a surprising link between homological algebra and Hadamard 
matrices jH] led to the study of cocyclic Hadamard matrices [17] . Hadamard matrices 
of many types are revealed to be (equivalent to) cocyclic matrices [18] . Among 
them, Sylvester Hadamard matrices, Williamson Hadamard matrices, Ito Hadamard 
matrices and Paley Hadamard matrices. Furthermore, the cocyclic construction is 
the most uniform construction technique for Hadamard matrices currently known, 
and cocyclic Hadamard matrices may consequently provide a uniform approach to 
the famous Hadamard conjecture. 
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The main advantages of the cocyclic framework concerning At by At Hadamard 
matrices may be summarized in the following facts: 

• The test to decide whether a cocyclic matrix is Hadamard runs in 0(t 2 ) time, 
better than the 0(t 3 ) algorithm for usual (not necessarily cocyclic) matrices. 

• The search space is reduced to the set of cocyclic matrices over a given group 
(that is, 2 s matrices, provided that a basis for cocycles over G consists of s 
generators), instead of the whole set of 2 16 *" matrices of order At with entries 
in {-1,1}. 

In pQ was shown that the cocyclic technique can certainly be extended to handle 
the maximal determinant problem at least when n = 2 (mod 4). More concretely, 
the study focused on cocyclic matrices over the dihedral group of 2t elements with t 
odd. Based on exhaustive and heuristic searches, three algorithms for constructing 
cocyclic matrices with large determinants were provided. 

In this paper we are interested in embedding (cocyclic) submatrices of orders 
2t with large determinants in certain cocyclic Hadamard matrices of orders At. If 
these determinants attain the largest possible value, we are embedding D-optimal 
designs. Also, we discuss the relation between the existence of these submatrices and 
the growth factor for these Hadamard matrices. 

In Section 2, an algebraic formalism (in terms of cocycles) to describe two com- 
binatorial operations on a matrix (eliminate and add certain rows and columns) is 
provided. As a consequence of this formalism a method arises for embedding (co- 
cyclic) submatrices of orders 2t with large determinants in certain cocyclic Hadamard 
matrices of orders At. In Section 3, we connect the existence of specific matrices em- 
bedded in cocyclic Hadamard matrices of order 20 with the values of the pivots that 
appear when we perform Gaussian elimination with complete pivoting on them. The 
last section is devoted to conclusions and future work. 

Notation. Throughout this paper we use — for —1 and 1 for +1. We write H for a 
Hadamard matrix and Dj for a D-optimal design of order j. The notation Dj £ H 
means Dj is embedded in H. 

2. Cocyclic D-optimal designs embedded in Cocyclic Hadamard ma- 
trices. Assume throughout that G = {gi = 1, g2, ■ ■ ■ ,g n } is a multiplicative group, 
not necessarily abelian. Functions ip: G x G — > (—1) = Z2 which satisfy 

(2.1) i>(9i,9j)i>(9i9j,9k) = i>(9j,9k)ip(9i,gj9k), ^9i,9j,9k G G 

are called (binary) cocycles (over G) [21] . A cocycle is a coboundary d<j> if it is derived 
from a set mapping G — s> (—1) by d<p(a, b) — 0(a)<^(6)0(ai>) _1 . 

A cocycle ip is naturally displayed as a cocyclic matrix (or G-matrix) that is, 
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the entry in the (i, j)th position of the cocyclic matrix is ip(gi,gj), for all 1 < i,j < n. 

A cocycle ip is normalized if ip(l,gj) = ip(gi,l) = 1 for all gi,gj G G. The 
cocyclic matrix coming from a normalized cocycle is called normalized as well. Each 
unnormalizcd cocycle ip determines a normalized one —ip, and vice versa. Therefore, 
we may reduce, without loss of generality, to the case of normalized cocycles. 

The set of cocycles forms an abelian group Z(G) under pointwise multiplication, 
and the coboundaries form a subgroup B(G). A basis B for cocycles over G consists 
of some elementary coboundaries di and some representative cocycles, so that every 
cocyclic matrix admits a unique representation as a Hadamard (pointwise) product 
M = Mq^ o . . . o Md iw o i?, in terms of some coboundary matrices Mg i . and a matrix 
R formed from representative cocycles. 

Recall that every elementary coboundary dd is constructed from the characteristic 
set map dd- G — > {— 1, 1} associated with an element gd G G, so that 



Remark 2.1. (J^ Lemma 1]) 

In particular, for d ^ 1, every row s ^ {1, d} in Mg d contains precisely two —Is, 
which are located at the positions (s,d) and (s,e), for g e = g^gd- Furthermore, the 
first row is always formed by Is, while the d-th row is formed all by —Is, excepting 
the positions (d, 1) and (d, d). 

Although the elementary coboundaries generate the set of all coboundaries, they 
might not be linearly independent (see [5] for details). 

Let G r (M) (resp. G C (M)) be the Gram matrix of the rows (resp. columns) of 



dd{gi,gj) = Sd(g l )Sd{gj)Sd(g l gj) for 5 d (gi) 




M, 



G r (M) = MM T , (resp. G C (M) = M T M). 



The Gram matrices of a cocyclic matrix can be calculated as follows. 



Proposition 2.2. (JTR, lemma 6.6]) 
Let be a cocyclic matrix, 
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If a cocyclic matrix My is Hadamard, we say that the cocycle involved, tp, is 
orthogonal and My is a cocyclic Hadamard matrix. The cocyclic Hadamard test 
asserts that a normalized cocyclic matrix is Hadamard if and only if every row sum 
(apart from the first) is zero |17j . In fact, this is a straightforward consequence of 
Proposition 12.21 



Analyzing this relation from a new perspective, one could think of normalized 
cocyclic matrices meeting the bound (|l.ip as normalized cocyclic matrices for which 
every row sum is zero. Could it be possible that such a relation translates somehow 
to the case n = 2 (mod 4)? We now prove that, in fact, the answer to this question 
is affirmative. 

A natural way to measure if the rows of a normalized cocyclic matrix M = [my] 
are close to sum zero, is to define an absolute row excess function RE, such that 



RE(M) = £ 



3=1 



This is a natural extension of the usual notion of excess of a Hadamard matrix, E(H), 
which consists in the summation of the entries of H . 

With this definition at hand, it is evident that a cocyclic matrix M is Hadamard 
if and only if RE{M) — 0. That is, a cocyclic matrix M meets (jl.ip if and only if 
RE{M) is minimum. This condition may be generalized to the case n = 2 (mod 4). 

For the remainder of the paper t denotes an odd positive integer. 

Proposition 2.3. fl]/ Let M be a normalized cocyclic matrix over G of order 
n = 2t. Then 
RE(M) >2t-2. 

But we may go even further. Having the minimum possible value 2t — 2 is a 
necessary condition for a cocyclic matrix M to meet the bound (|1.2I) . 



Proposition 2.4. FJ]/ If a cocyclic matrix M of order n — 2t meets the bound 
iTO) . then RE(M) = 2t-2. 

Unfortunately, although having minimum absolute row excess is a necessary and 
sufficient condition for meeting the bound (II. ip . it is just a necessary (but not suffi- 
cient, in general, see [1] Table 5] ) condition for meeting the bound (|1.2p . But there 
are some empirical evidences that matrices having minimum absolute row excess cor- 
respond with matrices having large determinants, see Table 2.1., page 11. 

From now on, we fix G = T>2m as the dihedral group with presentation (a, b: a m — 
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b 2 = (ab) 2 = 1), with ordering {l,a, ...,a m 1 ,b 7 ab, . . . ,a m 1 b} and indexed as 
{1, . . . , 2m}. A basis for cocycles over D 2m consists in (see [HE]): 

• Let nr be an odd positive integer. 

B= {3 2 ,...,<9 2 m-l,^2}- 

• Let m be an even positive integer. 

B = {d 2 , . . . ,<9 2m -2, ^1,^2,7}- 

Here di denotes the coboundary associated with the i t/l -element of the dihedral 
group T>2m, that is a l ~ x ( modm )&L-^-J _ And f3 1} /3 2 and 7 are the representative cocycles 
in cohomology, i.e. the cocyclic matrices coming from inflation are Mp ± = J m eg) 



1 1 

1 - 



and Mp 2 = 



1 1 

f - 



5 J m . We use A®B for denoting the usual Kronecker 



product of matrices, that is, the block matrix whose blocks are aijB. 

for the m x m 



The transgression cocyclic matrix M 7 is M 1 
matrices A m — iflij) and B m — (&y) where 



4 A 
B m B 

■11 



- 1 i + j > m + 1 J — 1 i < j 

1 otherwise ' y \ 1 otherwise ' 

It has been observed that cocyclic Hadamard matrices over the dihedral group mostly 
use M/3 2 oMj and do not use Mp 1 (see [3 [12]). In the sequel, we consider only cocyclic 
Hadamard matrices of this form Mq h o . . . o Mg im o Mp 2 o M 1 , 

In what follows, the goal is to provide an algebraic formalism (in terms of co- 
cycles) to describe two combinatorial operations on a matrix: the first consisting in 
eliminating and the second in adding certain rows and columns. 

Remark 2.5. T>2t is trivially embedded as a subgroup of T>4t, the dihedral group 
of it elements. Concretely, if V it = (a, b: a 2t = b 2 = (ab) 2 = 1) then V 2t = {a 2 ,b) C 

Proposition 2.6. Let be a cocyclic matrix overT)^, then the 2t by It matrix 
obtained by eliminating from the rows and columns indexed with an even number 
is a cocyclic matrix over 2? 2t and we denote it as . 

Proof. On the one hand, taking into account the ordering fixed above, 2?4 t = 
{1, a, . . . , a 2 *" 1 , 6, a6, . . . , a 2t_1 6}, the rows and columns in indexed with an odd 
number correspond with {1, a 2 , . . . , a 2t ~ 2 , b, a 2 b, . . . , a 2t ~ 2 b} = (a 2 , b). 

On the other hand, if *0| 2 b} denotes the restriction of i/> to the subgroup (a 2 , b), 
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then ip\ a satisfies (|2.1j) for G = (a 2 , b) = 2?2t since tp satisfies (|2.1[) for G — T>4 t ■ 
In other words, 2 b} i s a cocycle for G = {a 2 , b) and My — a . □ 

Lemma 2.7. Let M^ 1 and My 2 6e cocyclic matrices overT>^ t then 

w/iere M^.^ = o M^, 2 . 

Proof. It is a straightforward consequence of the pointwise multiplication. □ 

From now on, B = {<9 2 , • ■ • , 9 4t _ 2 , #2, 7} and B = {9 2 , • . • , fot-i,^} denote a 
basis of cocycles for T>^ t and 2?2t, respectively. 

Lemma 2.8. The following identities hold: 

{ Jit i even _ \ J 2 t i = 1 
Md ' = { M B^ iodd > = I * = 2 

and 

2 

i=l 

Proof. The identities above follow by direct inspection. □ 

Given My a 2?4 t -matrix. The following result describes, in terms of cocycles, the 
unique 2?2t-matrix obtained by eliminating from My the rows and columns indexed 
with an even number. 

Theorem 2.9. Given a V it -matrix 

= M? 2 0...0 M" 4t - 2 o M)} o M k r? o 

V 02 044—2 PI P2 7 

where (a.2, ■ ■ ■ > k\,k2, k^) denotes a concrete it-upla with entries or 1. Then 

= M? 2 0...0 M" 4t " 2 o M.' 1 o M k H 2 o M k3 

V 02 Oit-2 PI P2 7 

(t-1 
— 
11 9 2 ; d 2 t-2i+i 
i=l 

Proof. It follows from Lemmas 12.71 and 12.81 □ 




Embedding cocyclic D-optimal designs in cocyclic Hadamard matrices 



9 



Now, fixed a X^t-matrix. The whole set of 2?4t -matrices constructed by 
adding certain rows and columns to Mi is provided in the next theorem. 

Theorem 2.f0. Given a V 2 t -matrix 
Then any V '^-matrix of the form: 

\ k2 

= II M&oMfcoMfo n Mg u _, o Mg it _ u+1 o Mg^o. . .oMg^oM^, 

i=l \ i=1 J 

with on, k\ and k 2 taking values from {0, 1}, satisfies that = Mr. 

Proof. Using Theorem l2.9[ it is easy to check that — M^. □ 

Using our algebraic formalism we are able to give a method for embedding a 
2?2t-matrix in the rows and columns indexed with an odd number of a I?4 t -Hadamard 
matrix whenever it is possible. Although, theoretically this method provides a solu- 
tion, from the practical perspective it is only appropriate for numerical calculations 
in low orders, because in the worst case it needs to check 2 2t possibilities. Hence, 
properties providing some cut down in complexity and some heuristic are generally 
needed. 



Algorithm 2.11. Embedding 2?2*-niatrices in X>4 t -Hadamard matrices. 
Input: a 2?2t-niatrix. = Mg, o . . . o Mg. o M^. 

Output: a T>4 t - Hadamard matrix which contains embedded (that is, = 
M£), if such matrix exists. 

Step 1. Calculate = Mg 2ii _ 1 o . . . o Mg 2im _ 1 o Mp 2 . 
Step 2. Calculate all possible combinations 

\ k 

M,f, = \{ o M* o Yl Mg u _ 1 o M d4t _ it+1 o M 4 , 

i=l \i=l j 

where k and at may take the values or 1. 

Step 3. If there exists a (ai, . . . , ot-xt— i; k) such that RE(M c j ) ) = 0, then 
is Hadamard. Otherwise, such matrix does not exist for M^. 

Verification: By construction, Mq is a cocyclic Hadamard matrix over D^t and 
M$ = Mr (see Theorem HHflJ . 
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In the sequel, we describe an algorithm looking for 2?2*-matrices with large deter- 
minant embedded in a Hadamard matrix. Actually, the output is a Hadamard 
matrix where the £>2t-niatrix, M^, is embedded in a such a way that = . 
We give two strategies that are both based on exhaustive searches. The first one 
needs to calculate 2?2t-matrices with large determinant (by the approach given in [I]) 
and use Algorithm 12.111 for detecting if it is embedded in a Hadamard matrix. The 
second one needs to calculate 2?4t-Hadamard matrices, M$, (using [3]) and check the 
determinant of M^> ■ 

Algorithm 2.12. Search for 2?4 t -Hadamard matrices where a X^t-matrix with 
large determinant exists embedded in it. 

Input: an odd positive integer t. 

Output: a 2?4t-Hadamard matrix where ^ t ^2)(2t-2y- 1 — K ' ^ sucn matrix exists, 
with 0.85 <K<1. 

Approach 1 
0^0 

S <— The list of 2?2t-matrices satisfying that ( 4t _ 2 )(2t-2) t - 1 — K 
while S is not empty { 

1. Choose a matrix Mx in S. 

2. S^S\{M^}. 

3. Check (using Algorithm 12. lip whether exists a 2?4t-Hadamard matrix M^p 
such that M.0 = M^. If not, go to 1; otherwise Q -s— M^, 5 = 0. 

4. End while. 

} n 

Approach 2 

S <— The list of 2?4 t -Hadamard matrices 
while S is not empty { 

1. Choose a matrix in S. 

2. 5<-S\{Afy}. 

3. Check whether ^^^{^t-ly- 1 — K ' ^ n0 ^' §° ^° ^' otherwise Q <— M^, 
5 = 0. 

4. End while. 
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t 


#(M,/,) 


det(AL h )/2 2t - 1 




RE 


R 


3 


72 


5 


36(6) 


4 


1 




4 


36(9) 


4 


0.8 






144 


100(25) 


8 


1 


5 


1400 


125 
81 


1200(100) 
100(50) 


8 
8 


0.86 
0.56 






9477 


392(196) 


12 


1 






8405 


2352(294) 


12 


0.887 






7569 


2352(294) 


12 


0.799 


7 


11368 


4096 


392(196) 


12 


0.432 


2197 


1176(294) 


20 


0.232 






845 


1764(294) 


20 


0.0892 






841 


1176(294) 


20 


0.0887 






576 


1764(147) 


20 


0.061 



Table 2.1 

Determinant spectrum for T>2t -matrices embedded in Hadamard T>4t -matrices, t = 3,5, 7. 



} n 



For 3 < t < 9 odd, we have performed an exhaustive search looking for the pos- 
sible values of the determinant of £>2t-matrices embedded in 2?4 t -Hadamard matrices 
using Algorithm 12.121 with < k < 1, and have been displayed in Tables 2.1 and 
2.2. All the calculations have been worked out in Mathematica 4.0, running on a 
Pentium IV 2.400 Mhz DIMM DDR266 512 MB. 

In Tables 2.1 and 2.2, #(M^,) denotes the number of 2?4 t -Hadamard matrices 
(with representative cocycle Mp 2 o M 7 ). The second column of the table shows the 
different values for the determinant of M^, whereas the third displays about the 
frequency of £>4 t -Hadamard matrices with the same value for det{M^) and between 
parenthesis is the number of different 2?2t-matrices embedded, the fourth column 
gives the absolute row excess of M^. Finally, the last column informs about R = 

{it-2){2t-2) t - i ( wmcn is called efficiency of the design in |27|). 

Some interesting properties can be observed in these tables. First, for t = 3, 5, 7 
and 9 there is just one value of R above 0.9 (when R — 1 we got D-optimal designs) 
and some of them above 0.85. Second, for every pair of Hadamard matrices M J p 1 
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94 
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130248 


240448 


3888(486) 


24 


0.216 


9 


201977 


1944(486) 


24 


0.181 






198005 


3888(972) 


24 


0.178 






186624 


7776(972) 


24 


0.168 






179776 


1944(486) 


24 


0.161 






155236 


3888(1944) 


24 


0.139 






126025 


5832(972) 


24 


0.113 






90593 


972(486) 


32 


0.081 






88445 


1944(972) 


24 


0.079 






59049 


5832(162) 


32 


0.053 






20800 


5832(486) 


24 


0.019 






4096 


972(243) 


32 


0.004 









5832(486) 


32 






Table 2.2 

Determinant spectrum for T>\%-matrices embedded in Hadamard T>^Q-matrices. 



and M 4 , 2 with det{M^ 2 ) = det{M^) satisfy that RE{M^ t ) = RE(M^ 2 ). Third, 
the minimum RE corresponds with the largest values for the determinant for M^,. 
Taking into account this last property, a result about characterizing the 2?2t-matrices 
with minimum RE has interest. This characterization will be given here in terms of 
maximal n-paths, for this terminology we refer to [2]. 

Proposition 2.13. Let be a V 2t -matrix. The RE(MA = 2t - 2 iff 

1. decomposes as a combination of the form Mq o . . . o Mg o . 

2. The number of maximal s -paths for is either or for s = 2, .., t. 

Proof. The proof of this result follows just from the study of the distribution of 
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— 1 by rows in the elementary coboundaries (see [2]). Actually, considering a matrix 
N = AIq, o . . . o Mg, , Remark 12 . II implies that there is (necessarily) an even positive 
number 2/ s of —Is located at row s of N, 2 < s < 2t. Hence, 

RE(N) > At - 2. 

For t + 1 < s < 2t, taking into account the presentation of D^t 1 it may be readily 
checked that (a fc fe) _1 = a b. In this circumstance, the —Is entries located at row s 
are distributed in such a way that precisely f s of them occur through columns 1 to i, 
whereas the remaining f s occur through columns t + 1 to 2i. Furthermore, focusing on 
row s, any two coboundary matrices Mg. and Mg. either share their two —Is entries 
at row s (just in a unique row), or do not share any of them at all s. 

Consequently, attending to the form of , the summation of row s of Mr is 
iBMf = M §ti o...oM Stv> °Mg 2 . 

For 2 < s < t, attending to the form of and the fact that for every row s every 
maximal s-path determines two — l's entries in N (see [5])- Then, the summation of 
row s of ^ is 2 or —2 iff the number of maximal s-paths for M is either or £ti 
for s = 2, .., t. 

From above two equivalences, we conclude with the desired result. □ 

3. On the pivot structure of 2?20-Hadamard matrices. In 1988 Day and 
Peterson [7] proved that the growth factor of an n x n Sylvester Hadamard matrix is 
n. In [THl Example 6.2.4] it is shown that the Sylvester Hadamard matrix of order 
2" is cocylic over ZJ?. Taking into account the difficulty of proving Cryer's conjecture 
and the previous statements, the following question arises in a natural way: Is the 
growth factor of a £>4 t -Hadamard matrix its order? 

In this section, we will focus on 2? 2 o-Hadamard matrices and will extract some 
information on their pivot patterns from the D-optimal designs (and other specific 
sub-matrices with concrete determinants) that exist embedded in these matrices. 

It is well-known that there is a unique D-optimal design up to equivalence in each 
order up to n = 10 23 . 

In Table 2.1, the D-optimal design of order 10 appears (see !3.1|) embedded in 100 
different 2?20-Hadamard matrices. Using Magma V2.ll we found just one inequivalent 
Z^o-Hadamard matrix among the 100 (in higher orders, see [1] for an inequivalence 
criterion for X>4 t -Hadamard matrices). Based on exhaustive searches for minors of 
order 10 embedded in Hadamard matrices of order 20, we also identified the D- 
optimal design of order 10 in the other two inequivalent Hadamard matrices of order 
20. So, we can state the following result: 
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Lemma 3.1. The D-optimal design of order 10 is embedded in all Hadamard 
matrices of order 20. 

Corollary 3.2. For a Hadamard matrix of order 20, it holds that H(10) can 
take the value 144 • 2 9 . 



(3.1) 



1 



1 1 

- 1 
1 1 

- 1 

- 1 
1 1 
1 1 

- 1 
1 - 



1 
1 
1 

1 

-1 
-1 
1 



11- 



Proposition 3.3. Given H a Hadamard matrix of order 20, there exists an 
equivalent CP matrix H' where D\q appears located at the top left position in H' . 

Proof. Let H be an equivalent matrix to H such that Diq appears located at the 
top left position in H. Notice that Dio (see (|3.1|) is feasible. Our aim is to prove that 
D 10 (i)=H'(i), i = l,...,10. 

• Hi — 1,2,3,4,5,6 and 10 then the identity holds because -Dio(*) takes the 
maximum determinant of any (—1, l)-matrix of order i (see [25]). 

• It well-known that D 6 (£ D 7 ( [26l [24] ) , hence the maximum determinant of 
any (—1, l)-matrix, A, where Dq appears located at the top left position, 
is lesser than or equal to 512 (the second greatest, see the spectrum of the 
determinant function Due to D w (7) = 512 then D w (7) = H'(7). 

• We extend the 7x7 matrix -Dio(7 x 7) located at the top left position in Dio 
to all possible 8x8 matrices of the form 



where the elements * can be ±1. Computing the determinant of these 2 
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8x8 matrices, we have that it is always lesser than or equal to D\q{7) = 2560. 
Therefore, D w {8) = H'(8). 

• Taking into account that Dq ^ -Did then the maximum 9x9 minor of Diq 
is lesser than or equal to 12288 (the second greatest, see the spectrum of the 
determinant function [25]). Due to D 1Q (9) = 12288 then L>i (9) = H'(9). 

Now, taking H' as the CP extension of H, the proof follows. □ 

Corollary 3.4. The pivot pattern (1, 2, 2, 4, 3, 10/3, 16/5, 5, 24/5, 6) for the first 
ten pivots appears in the three classes for Hadamard matrices of orden 20. 

In a recent search that we have performed, we have found this result: 

• There is a CP matrix H equivalent to the following 2?20 _ H a damard matrix 
My = Mg^oMa^M^oM^oM^oM^oMftoM, such that #(10) = 125-2 9 . 

Although we didn't find any result in the literature asserting that if the existence 
of a submatrix with large determinant is proven for a matrix A, then we can indeed as- 
sume that it always appears in the upper left corner for some CP matrix A' equivalent 
to A. It seems to be true at least when this submatrix reaches the largest determinant 
(see [Ml p. 1763]). Also, the above result agree this for other large determinant (the 
second largest in Table 2.1). 

4. Conclusions and further work. In this paper we have described a method 
for embedding a 2?2t-matrix in the rows and columns indexed with an odd number of 
a I?44-Hadamard matrix whenever it is possible. If this 2?2t-matrix has a determinant 
attaining the largest possible value, we get a D-optimal design. This method relies 
on two combinatorial operations on a cocyclic matrix: eliminate and add certain 
rows and columns. The idea behind this approach has been to translate these two 
combinatorial operations into a pure algebraic framework (concretely, in terms of 
cocycles). Finally, our study has provided some information about the pivot values 
when Gaussian elimination with complete pivoting is performed on Z?2(rHadamard 
matrices. 

Our next goals are: 

• Study the relationship between the RE and the values of the determinant for 
2?2t-matrices. 

• Design heuristic searches based on RE for Algorithm 12. 121 

• Study if My satisfies that ^^(^t-ly- 1 — ^-85 implies that a CP matrix M 
equivalent to My exists such that M(2t) = det(My). 

• Study the pivot structure of 2? 4t -Hadamard matrices. 

• Design an "efficient" method to construct 2?4 t -Hadamard matrices from T>2t~ 
matrices. 
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